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Recently Gouy rotation was observed with focused non-relativistic electron vortex beams. If the 
electrons in vortex beams are very fast we have to take into account relativistic effects to completely 
describe the Gouy phase on them. Exact Hermite-Gaussian solutions to the Klein-Gordon equation 
for particle beams are obtained here that make explicit the 4-position of the focal point of the beam. 
These are Bateman-Hillion solutions with modified phase factors to take into account the rest mass 
of the particles. They enable a relativistic expression for the Gouy phase to be determined. It is in 
fact shown all the solutions are form invariant under Lorentz transformations. It is further shown for 
the exact solutions to correspond to those of the Schrodinger equation the relative time between the 
focal point and any point in the beam must be constrained to be a specific function of the relative 
spatial coordinates. 
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Introduction - The Gouy phase of matter waves in 
the non-relativistic context have been explored firstly in 
00, followed by experimental realizations with Bose- 
Einstein condensates Q, electron vortex beams [ 4 ] and 
astigmatic electron matter waves using in-line holography 
0. Also, it was showed that to improve the accuracy in 
imaging the dynamics of free-electron Landau states the 
diffractive Gouy phase rotation has to be reduced Q. 


Electron vortex beams recently have been theoretically 
and experimentally explored since they can be applied to 
improved electron microscopy of magnetic and biological 
specimens 0 0 . Production of electron vortex beams 
with high-energy and orbital quantum number seems fea¬ 
sible [9j. The Angular Momentum and Spin-Orbit inter¬ 
action in the relativistic electron vortex beams was stud¬ 
ied in terms of the Dirac equation in 0. The Gouy 
phase for this case, which was not treated in 0 , is com¬ 
pletely explained only if relativistic effects be considered 
and in this letter we explore these possible effects in the 
Gouy phase. The experimentally observed solutions for 
electron vortex beams with carry out Orbital Angular 
Momentum (OAM) could be expressed in terms of the 
solutions that will be obtained here. In order to avoid 
additional complications we do not take into account the 
spin of particle and we will solve here the Klein-Gordon 
equation. 


Bateman [11] discovered a class of exact solutions 
to the linear wave equation over a hundred years ago. 
Hillion [Hill later complexified these solutions for ap¬ 
plication to wave packet and wave beam problems. We 
present exact Bateman-Hillion solutions to the Klein- 
Gordon equation for the wave function ’I' ( x in Hermite- 
Gaussian beams where X/j, denotes position in Minkowski 
space and p, = 0,1, 2,3. Here we characterize the quan¬ 
tum dynamics of a relativistic particle enabling us to bet¬ 
ter understand the quantum behavior of electrons in high 
velocities, for example in relativistic electron vortex. As 


a consequence of the relativistic effects the Gouy phase 
for the matter waves has a different expression from the 
non-relativistic counterpart that must be taken into ac¬ 
count, especially, when we focus on a relativistic particle 
beam. Also, we present a method for using these solu¬ 
tions to calculate the properties of beams for continuous 
wave sources. 


The origin of the coordinate system for a beam is usu¬ 
ally at the center of the focal point. It follows the spatial 
coordinates x\{i = 1,2.3) of this point are zeroed out 
but translating the beam will make them explicit. For 
relativistic solutions, it must be recognized xI is part of 
a 4-vector such that the 4-potential is also dependent on 
the time coordinate Xq. 

The fact 4/ depends on two 4-position vectors creates 
a problem familiar from the treatment of two interacting 
relativistic particles 00 that the field cannot evolve in 
two independent time coordinates. The known solution 
to be applied here is to use a Dirac delta function S [/(G)] 
to impose a relationship /(G) = 0 between the relative 
coordinates G = x u — x^. In classical electrodynamics, 
this idea leads to the derivation of the Lienard-Wiechert 
potentials 0 for the field experienced at one point owing 
to the presence of a point charge at another. The delta 
function notation for eliminating the relative time is also 
used in quantum electrodynamics 171. 


For non-relativistic beams the wave function T can also 
be calculated from the Schrodinger equation 0 . Solu¬ 
tions to the Schrodinger equation are harmonic in time 
but 14/1 is independent of time. It is usual therefore to 
normalize 41 in a 3-dimensional constraint space rather 
than over all Minkowski space. For the exact solutions to 
be developed in this letter |4/| does depend on time. The 
intent though is still to normalize them in a 3-dimensional 
constraint space using the 5 [/(G)] condition to define it. 

Hermite-Gaussian Beams - Consider a beam of parti¬ 
cles each having a rest mass mo, a 4-position x M = ( Xi , t ) 
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and a 4-momentum p^ = (pi,E). The Klein-Gordon 
equation for the wave function ^(a^) representing each 
of the particles in Minkowski space can be expressed as 


where 


{E 2 + c 2 p 2 )’F = tooc 2 4/ 

( 1 ) 

„ id - id 

( 2 ) 

h dxi ’ e hdt 


are the momentum and energy operators, h is Planck’s 
constant divided by 2-7T and c is the velocity of light. 

Here we will find an exact solutions 'k mn to the Klein- 
Gordon equation m for Hermite-Gaussian beams moving 
in the 23 -direction. It will be assumed this takes the 
Bateman inspired form 


^ mn *&mn (Cl) £ 2,£3 + ct) exp[i(fc 3 2 3 - cot)\, (3) 


where 


M&, r) = w 0 ^l+ ) 3 (11) 

is the beam radius, wq = w( 0 , 0 ) is the beam waist and 

ffmn(C 3 , t ) = (l+m + n) arctan ^ , ( 12 ) 

is the Gouy phase of a relativistic quantum particle. We 
observe that it depends of the longitudinal direction 23 
as the non-relativistic case but also it depends of the time 
and the speed of light. The Gouy phase is a necessary 
feature of our solutions that are relativistic and can be 
normalized. We observe that the relativistic effect in the 
Gouy phase can be negligible only if the evolution occur 
in a scale less than 1 ns and the longitudinal distance is 
large than 1 m. 

It is proposed the normalizing constant C mn in 
can be determined from the condition 


£i = 2 i- 2 f, T = t-t r 


(4) 


gives the position of each point 2 M relative to the cen¬ 
ter of the focal point of the beam 2 )), k 3 is the wave 
vector component along the direction of the beam, u> is 
the angular frequency and are scalar functions. The 
positive integers m and n indicate the mode of the beam. 

Inserting eq. ([3]) into the Klein-Gordon equation dTJ) 
gives 


' d 2 d 2 
dx\ + dx\ 



u_d\ 
c 2 dt) 


*&mn 


= 0 , 


(5) 


having spotted 


cf_ 

dx 2 


1 3 2 \ 

Jdt 2 ) 


4 *mn 


= 0 , 


( 6 ) 



\^mn\ 2 S[f(^)}d^=l, 


(13) 


where d 4 £ = d^id&d&dr and /(£ M ) is to be determined. 
For physical interpretation, it is clear a physical system 
cannot evolve in the two independent times t and r. The 
intent of the /(£ p ) = 0 constraint condition is therefore to 
eliminate the relative time in terms of the relative spatial 
coordinates. 

Lorentz Transformations - The wave function fl/ mn ( 2 M ) 
has been determined in eq. j 8 ]) to be an exact solution 
of the Klein-Gordon eq. © for a particle in a beam. 
The task ahead is to confirm these solutions are form 
preserving under Lorentz transformations. 

A 4-vector qp is relativistic if it preserves its form under 
the Lorentz transformation equations: 


owing to the dependence of $ mn on £3 and r in the linear 
combination £3 + ct. 

Equation J5]) is analogous to the paraxial wave equa¬ 
tion 


dx\ dx 2 3 dx 3 


(7) 


where the role of £3 in <I> p (£i, £ 2 , £ 3 ) is replaced by £3 + cr 
in $ (£ 1 , £2, £3 +ct) . It can therefore be solved analogously 
flit to give 


(5 _ CmnWO „ 

^ m.n. — -LJ-rt 


\/ 2 £i 


H„ 


V2& 


w 


exp 


w 


w 


+ ^)(£i +£|) 

£3 + ct — 2 ib 


l 9n 


where H m and H n are Hermite polynomials, 


b = 


1 

2 



w 


2 

03 


( 8 ) 

(9) 


/ , Vi 

q l = 9i + 7 — 

c 


7 Vjdj 
+ 7 c 


- do 


(14) 


, / V jdj \ 

do = 7(90 - —) 


(15) 


(see ref. [16j]) where Vi is the relative velocity between 
any two inertial reference frames and 7 = (1 — u 2 ) -1 / 2 . 
For current purposes qp belongs to the set of position 2 M , 
relative position £ M , wave vector k^ and 4-momentum p ll . 

To confirm the Hermite-Gaussian beam solutions are 
fully relativistic it will be necessary to investigate the 
Lorentz transformation of number of different quantities 
including the phase factor, the Gaussian function, the 
Hermite functions, and Gouy phase. For brevity, the 
analysis will be limited to the case of an observer that is 
moving parallel to the axis of the beam. 

The product fc M 2 M is Lorentz covariant implying 




( 10 ) 


(16) 
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The numerator and denominator in the Gaussian com¬ 
ponent of eq. (0 can be transformed separately to give 



(17) 


and 

£3 + CT> ~ 2lb ' = \l (6 + ct- 2 ib), (18) 
\ c + v 

having determined 


For paraxial beams, the solution may be assumed 
to take the form 


^mn = $ mn(6,6,T) exp 


.h 


{P 3 X 3 - E s t) 


where 


(27) 


E s 


_P L 

2 mo ’ 


(28) 


is non-relativistic kinetic energy. Inserting eq. m into 
eq. (l26l) gives 


6 ' = 



(19) 


d 2 ® S mn , 

dx\ dx\ h dt 


(29) 


from eq. m assuming the radius wq is a Lorentz invari¬ 
ant scalar. Putting these results together leads to 


(*3 + t) (£f+& 2 ) = (fc 3 + f) fe 2 + C|) 

£3 + ct ' — 2 ib ' £3 + ct — 2 ib 


( 20 ) 


Eqs. (fT71) through CU1) can now be used to transform eq. 
0 into ®' mn is equal to 


CmnW 0 jj ( a/ 2^1 


-H, 




exp 


(^+^)( 6 2 + 6 2 ) , 
£3 + ct' — 2 ib' 


^1 X 

(21) 

w' J 


/ 

^ 9 mn 1 

(22) 


where 


w 


/ 


w 0 \ 1 + 



2 


w 0 



f &+CT 

V 2 b 


2 


(23) 


9 mn 

(1 + m + n) 


= arctan 


(^) 


= arctan 


( & 5 2 )- 

(24) 


and C mn is a Lorentz invariant scalar. 

It follows from taking stock of all of these results that 
eq. 0 can be rewritten as 


^'mn = £ 2 ,6 + CT ') exp(tk’ 3 x' 3 ~ Uj't'). ( 25 ) 

On comparing eqs. © and mb it is concluded 'F mn is 
form invariant under Lorentz transformations and there¬ 
fore fully relativistic. 

Correspondence to the Schrodinger Equation - The ob¬ 
jective here is to determine the form of the constraint 
condition /(6) = 0 in eq. (1131) from the principle that 
exact solutions of the Klein-Gordon equation should cor¬ 
respond to solutions of the Schrodinger equation in the 
non-relativistic limit. The Schrodinger equation takes 
the form 


Equation (l29l) is analogous to eq. © where the role of r 
in $ s (£i,£2,t) replaces £3 + ct in $(6,6.6 +cr). 

The solution to eq. (ESI is 


S CmnWQ I \/ 2 £i 


(T)° — 

ran 


H n 


exp 


Hn 


V 26 


n „s 

,J 9 mn 


moWQ + 2ihr 


where 


w s (t) = Wq y 1 + 

is the beam radius and 


Y 

Wo *0 / 


x (30) 

(31) 

(32) 


gmuir) = (1 + m + n) arctan ( -^ ] , (33) 

\m 0 w£J 


is the non-relativistic Gouy phase for matter waves. 

It can now be shown that the Schrodinger and exact 
Klein-Gordon beam solutions are related to each other 
through the expression 


$ 


s 

mn 


$mn< 5(6 - U 3 T)dt, 


(34) 


in the non-relativistic limit. This result is most easily 
demonstrated for the case of a Gaussian beam 


$00 = 


C 00 b 


exp 


b + * 5(6 + cr) 

On setting, 

m 0 W37 

k ‘ = -r- 

eq. EH evaluates to 

Goo m 0 wl 


*(*3 + f)fo 2 +& 2 ) 


6 + ct — i2b 


7TtQC 2 7 


^00 — 


rriQWQ + 2 ihr 


exp 


-m 0 (6 +^) n 


m 0 Wg 


(35) 


(36) 


(37) 


mn I w mn 

dx\ dx\ 


d 2 ^mn , ^ m d ^mn 
dx\ h dt 


= 0 . 


(26) 


having assumed 7 ~ 1. The calculation for all the higher 
order modes follows similarly. 
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It has been shown that in order for exact Bateman- 
Hillion solutions of the Klcin-Gordon equation for 
Hermite-Gaussian beams of matter waves to correspond 
to solutions of the Schrodinger equation in the non- 
relativistic limit, they must be restricted to the constraint 
space 

/&*) = 6 - u 3 t = 0. (38) 

This result is straightforward to interpret for particles 
large enough to have classical properties since the expres¬ 
sion 6 — u 3 t = 0 describes the trajectory of a classical 
free particle along the axis of the beam. It is unexpected 
that smaller quantum particles also appear to be confined 
to this constraint space. 

Discussion - Exact solutions have been derived to 
the Klein-Gordon equation for matter waves Hermite- 
Gaussian beams using Bateman-Hillion functions. The 
solutions have been shown to preserve their form under 
Lorentz transformations. These exact solutions are not 
time-harmonic but it has been found they can be used to 
produce results comparable to time-harmonic solutions 
providing they are normalized and interpreted in a 3- 
dimensional constraint space. The guiding principle is 
time-harmonic and Bateman-Hillion solutions can only 
be compared if the calculations are performed in spaces 
that have the same number of dimensions. 

It has been shown that for Klein-Gordon and 
Schrodinger beams to be in correspondence in the non- 
relativistic limit, the Klein-Gordon solution must be re¬ 
stricted to the <5(6 — M 3 r) constraint space. The con¬ 
straint equation 6 — u 3 t = 0 describes the trajectory of 
a classical particle. In this sense, it is intuitive that par¬ 
ticles large enough to exhibit classical properties would 
belong to this constraint space but surprising that smaller 
quantum particles also appear to be restricted to it. 
When the Gouy phase of matter waves was obtained in 
[l], the constraint equation 6 — u 3 t = 0 was used to 
transform the Schrodinger equation analogous to paraxial 
equation for which the propagation direction of a quan¬ 
tum particle is said to be classical. As the energy asso¬ 
ciated with the momentum of the particles in the prop¬ 
agation direction is very high to include the relativistic 
effects we can consider, as in [lj, a classical movement in 
this direction which makes the Gouy phase be a function 
of propagation direction and the relativistic effect appear 
in the time coordinate, see eq. (fl2l) . If we do not consider 
that the momentum of the particles in the propagation 
direction is high we have the non-relativistic Schrodinger 
equation that is different from the paraxial equation but 
still the Gouy phase is present in the solutions of a free 
particle as a property of the wave behavior. However, it 
is a time dependent function instead of a function depen¬ 
dent of the propagation direction. 

In conclusion, we have shown that the Gouy phase of 


relativistic quantum particle can be obtained by trans¬ 
forming the Klein-Gordon equation to an analogue of 
paraxial equation which provides normalized and Lorentz 
invariant solutions. The relativistic Hermite-Gaussian 
solutions contain Orbital Angular Momentum (OAM) 
and form a complete basis. These results can therefore 
be applied to treat relativistic electron vortex beams. It 
is anticipated though for strongly converging or diverg¬ 
ing beams that <5(6 — u 3 t) would need to be replaced to 
a form such as <5(|6 — \ui\r ) for better correspondence 
to the trajectories of classical particles in the case Too is 
large. 


[1] I. G. da Paz, M. C. Nemes , S. Padua, C. H. Monken, and 
J.G. Peixoto de Faria, Phys. Lett. A 374, 1660 (2010). 

[2] I. G. da Paz, P. L. Saldanha, M. C. Nemes, and J. 
G. Peixoto de Faria, New Journal of Phys. 13, 125005 
( 2011 ). 

[3] A. Hansen, J. T. Schultz, and N. P. Bigelow, Confer¬ 
ence on Coherence and Quantum Optics Rochester (New 
York, United States, 2013); J. T. Schultz, A. Hansen, and 
N. P. Bigelow, Opt. Lett. 39, 4271 (2014). 

[4] G. Guzzinati, P. Schattschneider, K. Y. Bliokh, F. Nori, 
and J. Verbeeck, Phys. Rev. Lett. 110, 093601 (2013). 

[5] T. C. Petersen, D. M. Paganin, M. Weyland, T. P. Sim¬ 
ula, S. A. Eastwood, and M. J. Morgan, Phys. Rev. A 
88, 043803 (2013). 

[6] P. Schattschneider, T. Schachinger, M. Stoger-Pollach, S. 
Loffler, A. Steiger-Thirsfeld, K. Y. Bliokh, and F. Nori, 
Nature Communications 5, 1 (2014). 

[7] J. B. McMorran, A. Agrawal, I. M. Anderson, A. A. Herz- 
ing, H. J. Lezec, J. J. McClellan, and J. Unguris, Science 
331, 192 (2011). 

[8] J. Verbeeck, H. Tian, and P. Schattschneider, Nature 
(London) 467, 301 (2010); P. Schattschneider and J. Ver¬ 
beeck, Ultramicroscopy 111, 1461 (2011); J. Verbeeck, P. 
Schattschneider, S. Lazar, M. Stoger-Polach, S. Loffler, 
A. Steiger-Thirsfeld, and G. Van Tendeloo, Appl. Phys. 
Lett. 99, 203109 (2011); S. Lloyd, M. Babiker, and J. 
Yuan, Phys. Rev. Lett. 108, 074802 (2012). 

[9] I. P. Ivanov, Phys. Rev. D 83, 093001 (2011). 

[10] K. Y. Bliokh, M. R. Dennis, and F. Nori, Phys. Rev. 
Lett. 107, 174802 (2011). 

[11] H. Bateman, Proc. London Math. Soc. 8, 223 (1910). 

[12] P. Hillion, J. Math. Phys. 33, 2749 (1992). 

[13] A. P. Kiselev, A. B. Plachenov, and P. Chamorro-Posada, 
Phys. Rev. A 85, 043835 (2012). 

[14] A. Komar, Phys. Rev. D 18, 1887 (1978). 

[15] H. W. Crater and P. Van Alstine, Phys. Rev. D 36, 3007 
(1987). 

[16] M. Saleem and M. Rafique, Special Relativity Applica¬ 
tions to Particle Physics and the Classical Theory of 
Fields (Ellis Horwood, 1992). 

[17] R.P. Feynman, Quantum Electrodynamics (Addison Wes¬ 
ley, 1962). 

[18] A. E. Siegman, Lasers (University Science Books, 1986). 



